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Abstract 

The correlators of light-quark currents contain mass-singularities of the form \og{rn? / Q"^) . 
It has been known for quite some time that these mass-logarithms can be absorbed into 
the vacuum expectation values of other operators of appropriate dimension, provided 
that schemes without normal-ordering are used. We discuss in detail this procedure 
for the case of the mass logarithms m^log(m^/Q^), including also the mixing with 
the other dimension-4 operators to two-loop order. As an application we present an 
improved QCD sum rule determination of the strange-quark mass. We obtain ifis{^ 
GeV) = 171 ± 15 MeV. 
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1 Introduction 



The method of QCD sum rules, first introduced in has become a popular and powerful 
technique to study QCD in the low-energy, non-perturbative region. The starting point is 
the Operator Product Expansion (OPE) of current correlators at short distances, suitably 
modified to incorporate non-perturbative effects. The latter are parametrized by a set of 
vacuum expectation values of the quark and gluon fields entering the QCD Lagrangian. 
These vacuum condensates induce power corrections to asymptotic freedom, and are respon- 
sible for the rich resonance structure observed at low energies. The basic assumption here 
is the factorization of short and long distance effects. The former are associated with the 
Wilson coefficients in the OPE, and the latter with the vacuum condensates. While the 
Wilson coefficients are calculable in perturbation theory, to any desired order in the strong 
coupling constant, the vacuum condensates cannot be calculated analytically from first prin- 
ciples (this would be tantamount to solving QCD exactly). Instead, they can be estimated 
in the framework of lattice QCD, or extracted from experimental data in certain channels 
by means of the QCD sum rules themselves. Next, making use of the analyticity properties 
of the relevant Green functions, and invoking the notion of QCD-hadron duality, one relates 
the fundamental QCD parameters entering the OPE with a dispersive integral involving the 
hadronic spectral function. In this fashion, a relation between hadronic and QCD parameters 
is achieved. 

An important problem which must be addressed in connection with the factorization of 
short and long distance effects in the OPE is the appearance of mass singularities in the 
coefficient functions. They are actually a long-distance effect and thus their presence in the 
coefficient functions spoils the desired factorization. It is possible, however, to shift them 
into the vacuum condensates, provided one is willing to accept the existence of perturba- 
tive vacuum expectation values of operators. This is equivalent to giving up the customary 
normal-ordering prescription which, by definition, sets such contributions to zero. A detailed 
discussion of how this can be achieved is presented in Section 2, including the renormaliza- 
tion group improvement. As a phenomenological application of these results, we address in 
Section 3 the problem of determining the value of the strange-quark mass. Some time ago 
a redetermination of the strange-quark mass was performed in the framework of QCD sum 
rules, exploiting new developments in the theoretical and experimental understand- 
ing of the two-point function involving the strangeness-changing vector current divergence. 
This constituted an improvement over earlier determinations of 0-0. Of particular im- 
portance was the removal of logarithmic quark-mass singularities in two-loop quark-mass 
corrections of order 0{m'^) and 0{m'j^), achieved in [Q. However, a feature of some concern 
was the appearance of parametrically enhanced terms of order 0{l/as)- In this paper we 
remedy this problem, showing how these terms can be effectively avoided by employing a 
scheme without normal ordering. Working at the three-loop level in perturbative QCD, and 
including two-loop radiative corrections to the condensates, we obtain a new expression for 
the current correlator. This is then used in order to obtain an improved value of the strange- 
quark mass. 
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2 The Operator Product Expansion 



We will be concerned in the following with the vacuum expectation value of the following 
time ordered product 

T{q) = i J dxe^''^r(J(x)J^(0)) , (1) 

where J = daS'y'^u = mssiu, and the up- and down-quark masses are neglected. Except for a 
sign change, which will be given explicitly, all our results will hold also for the divergence of 
the strangeness-changing axial vector current J = daSj^^j^u = msSi'j^u. When sandwiched 
between vacuum states, the T-product (1) becomes the corresponding 2-point correlator 

V'(g',a„m,/i) = (0|r(g)|0), (2) 

where n is the renormalization scale. Note that the polarization operator is not renormalization- 
group (RG) invariant as the function {0\T{J{x)J^{0))\0) contains non-integrable singularities 
in the vicinity of the point x = 0. These cannot be removed by the quark mass and cou- 
pling constant rcnormalizations alone, but must be subtracted independently. We write the 
correlator (2) as an expansion in powers of 1/Q as 



HQ^) = ml{Uo{L,as)Q'' + mlU2{L,as)±^^^^^^^{muUu) 







+ E^^^(o.-) + W-')}, (3) 

where L — ln(/x^/(5^), and Q'^ — —q^. The upper sign in front of the coefficient function Cu 
corresponds to the scalar case and the lower one to the pseudoscalar case. The operators Oi 
are 

(y 

0,^^Gl,G-^\ O^^msss, 03 = m^ (4) 

TT 

and the explicit expressions of the functions IIo and 112 will be given later. The only terms 
in %Ij which arc not RG-invariant are the coefficient functions m^IIo and m^H2- They satisfy 
the nonhomogeneous RG equations 

/^^("^^Ho) = m^7g , ^^[m^Jl^) = mj-fm , (5) 

with 

= (7) 

The anomalous dimension 7^ has been given to three-loop order because the corresponding 
correction to IIo is a priori not negligible, and it will be taken into account in the next 
Section. The absorptive part of iplQ'^), being an observable quantity, is invariant under the 
RG transformations. Without any loss of generality we will work with the second derivative 
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ip"{Q^) = dPipiQ"^) / diQ'^Y , which can be seen from (3) and (5) to satisfy an homogeneous 
RG equation 

i^^r{Q') = o. (8) 

The high energy behavior of tp'^Q"^) in the deep euchdean region may be rehably evaluated 
in QCD by employing the operator product expansion, i.e. 



Q V(Q^ «s, fiis, fx) ] ■ Kq{Q'^, as, TUs, /i) 1 

2 

+ E77^ E i^.(Q^«„m„/.)(0|O,(^)|0) 

n W ) dimOi=n 



(9) 



We have explicitly separated the contribution of the unit operator from that of the operators 
with a non trivial dependence on the field variables. The coefficient functions (CF) Kq and 
Ki depend upon the details of the renormalization prescription for the composite operators 
Oi. The usual procedure of normal ordering for the composite operators appearing on the 
r.h.s. of the OPE (9) becomes physically unacceptable if quark mass corrections are to be 
included. This is already obvious for the unit operator, representing the usual perturbative 
contributions if normal ordering is used, because in general it contains mass and momentum 
logarithms of the form 



(ln^)"(ln 



Q 



(10) 



with n, a and (3 being non-negative integers. More specifically, one may write 



<^(Q^«„m„/i) 



n>0,Z>0 



111 



where L = ln(/i^/Q^), M = ln(/i^/m^), and the overscript NO is a reminder of the normal 
ordering prescription being used. The function Fii{L,M) corresponds to the contribution of 
the /-loop diagrams, and is a polynomial of degree not higher than /, in both L and M. Now 
it is obvious that one may not choose the normalization scale fi in such a way that for Q ^ m 
both M and L would be small. The mass logarithms signal that even in the framework of 
perturbation theory there are effects coming from large distances of order 1/m. Fortunately, 
it has been realized long ago 0,0, 0-0 that all the mass logarithms may be neatly shifted 
into the vacuum expectation values (VEV) of non-trivial composite operators appearing on 
the r.h.s. of (9), provided the latter are minimally subtracted. 

To give a simple example, let us consider the correlator (2) in the lowest order one-loop 
approximation. First, we use the normal ordering prescription for the composite operators 
which appear into an OPE of the time ordered product in (1). To determine the coefficients 
of the various operators, one possible method is to sandwich both sides of the OPE between 
appropriate external states. By choosing them to be the vacuum, only the unit operator 1 
will contribute on the r.h.s., if the normal ordering prescription is used. This means that the 
bare loop of Fig.l contributes entirely to the coefficient Kq in (9). A simple calculation gives 
pi]| (in the sequel we neglect all terms of order and higher ): 



<^(g^ m„ /x) = Q''riQ'\ rn,, a,, /i) |„ 



Q 



2m^ 



M) 



■ (12) 
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This coefficient function contains mass-singularities (the M-term) . On the other hand, if one 
does not follow the normal ordering prescription, then the operator rrisSS develops a non- 
trivial vacuum expectation value even if the quark gluon interaction is turned off by setting 
as = 0. Indeed, after minimally removing its pole singularity, the one loop diagram of Fig. 
2 leads to the following result |0 



(0|.-.|0r=-^4iviln4 + l) . (13) 

By inserting this into (9), the new coefficient function Kq can be extracted, with the result 

3m? 



Kn 



(14) 



The mass logarithms are now completely transferred from the CF Kq to the VEV of the 
quark operator (13)! The same phenomenon continues to hold even after the as corrections 
are taken into account for (pseudo) scalar and (pseudo) vector correlators, independently of 
their flavour structure |n|]. The coefficient functions of the non-trivial operators will also 
depend on whether or not normal ordering is employed. 



The underlying reason for this was ffist established in |T^. Here it was discovered that if 
the minimal subtraction procedure is scrupulously observed [|13| then no CF may depend on 
mass logarithms in every order of perturbation theory, irrespectively of the specific model, 
and/or OPE at hand. This implies that all the mass logarithms log M in (11) go over into 
the "condensates", where they are hidden among various non-perturbative contributions. 
This remarkable property leads to the possibility of using the standard RG techniques to 
study mass effects in the framework of QCD sum rules without interference from unwanted 
mass singularities. On the other hand, the above nice features of minimal subtraction come 
at a price: when schemes without normal ordering are employed, then the renormalization 
properties of composite operators and CF's become more involved. This may already be 
observed in our one-loop example. Indeed, as a consequence of (8) one can immediately infer 
that the CF K^^ is RG invariant and hence, should obey the equation 

f^^Krm = o. (15) 

This equation is satisfied trivially for (12) but not for Kq as expressed by (14)! The reason is 
that the operator rrisSS ceases to be RG invariant in the world without normal ordering. The 
vacuum diagram of Fig. 2 has a divergent part which has to be removed by a new counterterm 
proportional to the operator m^l. In other words, rrisSS begins to mix with the "operator" 
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To lowest order, the corresponding anomalous dimension matrix reads 

(16) 





The nonvanishing off-diagonal matrix element describes the mixing of the two operators 
under renormalization and was obtained from the divergent part of the vacuum diagram in 



5 



Fig. 2. The diagonal matrix elements are just the anomalous dimensions of the respective 
operators in the usual normal-ordering scheme. The lower one is equal to —A'~f{as), where 
7(as) is the strange quark mass anomalous dimension which defines its running according 
to 



where 14 



d/x 

7(as) = 7i — + 72 ( 



+ 



with 7i = 2, 



72 

73 



101 5 

, 

12 18 ^' 



( 



1249 - 



2216 160^,^, 



n 



f 



(17) 
(18) 

(19) 
(20) 



and rif is the number of active light quarks. The running of the coupling constant as(/i) is 
determined by 



where 15 



+ 



(21) 
(22) 



with 



/92 



11 1 
51 19 

r + 12 

2857 5033 



1 Uf 

4 12 ^ 



325 

Uf 

18 ^ 54 



n}^ /32 . 



(23) 



The solutions of (17) and (21) can be written as 



asifi) 



2n 



'^^2^1nL^ 4 fPl 



(24) 



, , (, 271/32 2 / 71/^2 \ 1 

+^(-/5?73 + Phi + PIP212 + PiPhi - PiPsii - W1P21112 + PhD 

+ ^{PlP2i2 - P1P21112 + PhD + ^^i-PiPhi + PhD^j ■ (25) 

The expressions (24) and (25) are given to three-loop order for completeness. The simplified 
arguments so far will only make use of the corresponding one-loop results. Now one can see 
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that the operator rrisss acquires a scale dependence, which can be obtained from (16) and 
is given by 

mM^^) = mJsM + A £2 exp (-4 £^^^^ ^^M) <(/.o) 

A distinctive feature of this result is the appearance of inverse powers of 0. Note that 
in the approximation we have considered, the combination 

J, = (msSs)U) + —r^ — (27) 
^ 27r(47i + /5i) a,(/x) ^ ' 

is a RG invariant. It corresponds (but is not generally equal) to the RG- invariant combination 
rUsSS in the usual normal-ordering scheme. For simplicity, we will neglect for the moment the 
contributions of the dimension-4 operators and uu in the OPE (9). They will be added 
in later. The lowest-order coefficients of the operators rrisSS and at the scale = Q have 
the values: 

Cm.ss — 2 ' — Y^;^ ■ (28) 

We are now in a position to derive the RG improvement of the coefficient functions appearing 
in the OPE (9) when working in a scheme without normal ordering. To achieve this, one notes 
from Eq. (8) that the total contribution of the operators of dimension 4 is RG invariant and 
therefore we can choose freely the scale /i. Setting the renormalization scale fi = Q allows us 
to absorb all logarithms In //^/Q^ appearing in the CP Kq into the running coupling constant 
and the strange quark mass. On the other hand, the matrix elements of the operators at this 
scale can be expressed in terms of the same matrix elements at a lower scale /xq — 1 GeV^ 
with the help of the RG equation (16). Our first result for the RG improvement of the OPE 
(9), treated entirely within the minimal subtraction prescription, reads (for Uf = 3) 

Q-r{Q\ m„ ^^) ^ ^m^(Q) f 1 - 2^) 

+ ^((0|(-^-)(/^o)|0)--^ + -^j+0(l/Q). (29) 

This result has been essentially obtained for the first time in ^ (to two loop order). There 
are, however, a number of differences between its interpretation as given in (see also 



and [|TT|) and the point of view we will take in this paper. We comment briefly on these 
differences. In the vacuum expectation value of the RG invariant combination Ig in (27) 
was identifled with the (RG-invariant) product 771^(0 155^*^10) in the usual scheme using 
normal ordering. Thus, the vacuum matrix element (0|(msSs)(/i)|0) (non-normal ordered 
and RG-noninvariant) was represented, to the order we are working, as the sum of a RG- 
invariant part (of a nonperturbative origin, due to the spontaneous breaking of the chiral 
symmetry in QCD) and of a perturbative part (/i-dependent) which represents the sum of 



the leading mass singularities of the form ln^~^^ {fx^ / m^^ 



(0|(m,ss)(/i)|0) = (0|m,ss^«|0) - ^-^^ • (^O) 
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As it can be seen from (24) and (25) in the one-loop approximation, the second term van- 
ishes in the hmit fi —>■ oo and the distinction between the VEV of the operator rngSS in 
schemes with and without normal ordering disappears. It should be stressed that the above 
interpretation of the relation between normal ordered and non-ordered quark condensates 
relies heavily on an implicit assumption which is difficult to (dis)prove. Indeed, all purely 
perturbative contributions to (0|(msss)(/i)|0) were assumed to vanish in the limit /i — oo. 
Fortunately, even if the hypothesis fails it will only spoil the applicability of the scheme 
with normal ordering, but would have no effect on other renormalization schemes like e.g. 
the minimal subtraction prescription. The practical consequence of this approach is a large 
value of the mass correction of order m^, which is enhanced by the presence of one negative 
power of as{Q) (the second term in the second line of Eq.(29)). Note that in this approach 
there is no corresponding term containing l/as{fio) in Eq.(29), because it can be effectively 
combined with the contribution of the operator m^ss resulting in the RG-invariant VEV 
(0|Js|0). On the other hand, considering that a typical momentum transfer for QCD sum 
rules is of about 1 GeV, we will work with the quark and gluon condensates normalized at 
this "natural" scale fio = 1 GeV as our reference values. (For the case of the semihadronic 
decay rate of the tau lepton a similar approach has been suggested in [|16].) As mentioned 
above, this point of view is equivalent to the one taken in 0, provided the scale /xq would 
have been taken to infinity. Our choice of Hq somewhere around the characteristic momen- 
tum scale specific to the problem at hand (~ 1 GeV) helps to avoid the parametrically 
enhanced inverse powers of as- Indeed, as one can see from (24) to one-loop order, one has 
Ti /as{Q) — n/as{fio) = ^Pi ln(Q//io), which is not particularly large. 



3 Determination of the Strange Quark Mass 

We proceed now to include the contributions from the gluon operator and from the light 
quark condensate uu, working consistently to two-loop order. The coefficient functions Uq 
and 112 in (3) have, respectively, the three-loop and two-loop values (for an arbitrary 
renormalization scale /i, the use of the MS scheme is understood) 



no 



167r2 



+ 



TT 



12 - 

17645 
24 

,2 



TT 



353C(3) - 8 



511 3 
^/C(3) + —rif + -C(4) - 50C(5) 



-4n/C(3)ln^ + n,|ln^ + 117C(3)ln^ 



10801 , 



11 , 2 



2 



1061n^^ + nf-ln' 

g2 -5 



1 
3 



g2 2 g2 



(31) 



167r2 



TT 



-100 - 64 In ^ - 24 In^ ^ + 48C(3) 



(32) 



The three-loop result (31) was obtained in the works fl^ while the terms in (31,32) 
can be found in |T0, IT8[. Note that the first evaluation of Hq at three-loop level made in Ref. 
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|T9| proved to be erroneous ||20|. Unfortunately, that wrong result was then used in the work 
23| to find the light quark masses in the framework of the finite energy sum rules. 

The coefficient functions of the dimension-4 operators, evaluated at the scale fi = Q, are 

Cu = (1 + ^). (35) 
Cs = l|^(l + V^8C(3)-6)). (36) 

The leading order contributions to (33-35) were computed in The two- loop corrections 
to (33-36) were evaluated in the Refs. H, 0, H, i, pl and [|, ^ respectively. 

It is now a simple matter to derive the RG improvement of the no,2-terms in (3). Solving 
(5) with the boundary conditions (31)-(32) yields 

Kno)U = 

"4.(27, + A) «.(g) + + 4 + /^i)^ 

+ f + ^-8C(3)-2,.C(3))(^)^) (37) 
47r(27i + A) a,(/i) I ' vr H ^ V 



where 



ri=3-y + ^-^ = -2, ( = 3 ) (39) 
4 72 ^A /32 53 



9(7889 - 432C(3)) + 9(439 - 1824C(3))nj - 1195^2 
864(-57 + 2n/) 

5904C(3) -8011 



4896 



= 3) . (41) 



The RG improvement of the contribution of the dimension-4 operators in (3) requires the 
knowledge of their mixing under renormalization. The generalization of (16) to two-loop 
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order, by taking also into account the mixing with the gluon operator G^, reads 



nisSS 

4 

m:. 



/ 

V 



-a. 



das 













-470 
-47 



rrisSS 

\ 4 
V 



(42) 



Here /? and 7 were defined in (21) and (17) respectively, and 70 is the two-loop vacuum 
energy anomalous dimension, given by 



7o 



, 4 
1 + —a, 
Sit 



(43) 



The operator mgUU is RG invariant. The analogous anomalous dimension matrix which 
describes the mixing of the dimension-5 operators in schemes without normal ordering has 
been calculated recently, to one-loop order, in [^. We apply now the RG improvement of 
the contribution of the dimension-4 operators in (3). By taking advantage of the fact that 
their total contribution is RG-invariant, we choose /i = Q, where the coefficient functions 
are given by (33)- (36). The matrix elements of the operators Oi,2,3 can be scaled at /io — 1 
GeV with the help of (42), where they are known. This procedure leads to 



l^CjOj = (1 + — — — )(m,MM)o + - 



I 2 /3J vr 



+ 



/?2a.(/io) 



+ 



1 /II _ 7l 

2 I 6 2/?i 



Qs(Q) ^ 7i as(Aio) 



TT 



vr 



(02)1 



MO 



3 / as{Q). 47 _ 

47r(47i + A)««(g) \ TT ^^^'^12 4^ 



mt(Q) 



47r(47i + /3i)a,(/io) 



^h + ^ + -]l-.(/.o) 



«s(Q) 



47r2(47i + /?i)a,(/io; 



m^(/io) 11 - 



37/ 
A, 



(44) 



In order to keep the expressions within a reasonable size, we will replace here the various 
constants by their numerical values corresponding to = 3. At the same time, to help the 
reader who might want to reproduce our result, the nj-dependence which appears from other 
sources will be left explicit in the following. Thus, putting together (37) and (44) and taking 
two derivatives with respect to we obtain 



3m2(Q) llasiQ) /5071 



87r2Q2 



[1 



37r 
28a,(g)^ 
37r 



144 



f C(3)) (^' 
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+ 



ml{Q) 



4 TT 



2{msUu) ( 1 + 
1 



23a,(g)' 



I Mo 



16a,(/io) ^ 121a,(g) 



97r 



IStt 



4a,(/io) ^ 64a,(g)^ 



97r 97r 
TT 173' 



Ttt' 



' TT 155' 



64 a.(g) 4 



^as(/io) 72 / 2l7r2 ^^(yUo 



(45) 



A similar relation has been previously used in a QCD sum rule determination of the strange 
quark mass 0, where it was interpreted in the spirit of 0. As explained earlier, in the 
approach of ^ the normal-ordered strange quark condensate (times m^) is identified with 
the VEV of the RG-invariant combination Ig defined at one-loop level in Eq.(27). At two- loop 
level it has the form 



{msss){fi) + 



27r(47i a,(/i) 



77 



(46) 



Besides this, our result (45) differs from the one in p (see also |TT|) because there the mixing 
of the gluon condensate with the other operators of dimension 4 has been neglected. 



We perform the Borel transform L of ip'^Q"^), i.e. 



L[r{Q')] 



lim 



A simple calculation using the methods of gives 
3 m? 1 



X 



+ 



1 + 



91n^ 



^7r2 M2 [iln(MVA2)]-27i//9i 

/n^^4/32, , M2 4 

y-^^^(') + :^^"^"A^-A^ 

r2510167 1340^^3^^ 3^^,^^^ _ 



811n2(M2/A2) I 6561 9 



+ lnln(MVA2 

i2 



81 ^ 3 
1392641 



1109248 4352,,^,, , -.^ ^ . .'i, .2^ - ^ 
__ + _,)(l)J+,„^,„(MVA')^| 

4 



'M2 [iln(M2/A2)]-27i//3i 



1 + 



f-2..,(2) + 8|lnln^ 



/5i7i 



72 - 7i 



Pi 



2M6 [iln(M2/A2)]-27i//3i 



X 



A2 



5(^o) 



02 M2 4 



72 - 7i 



(47) 



A(/io 
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J 1 ml \ 155 /3i ^_^/67i \ 

77r2 2M6 [iln(MVA2)]-67i//3i I 24 2 A2 2 \ (3i j"^^' 

where 

^(/xo) = 2{msUU}o + -{—G }o{l + 

4 TT y y TT J 



N 46 . _ , 121 2s 64 . _ . 
= -:r{'msUu}o + —{—G }o + —{msSs}o 
o /^z TT y 

4/ ^ 64 TT 519 a,(/io) ^ 
2l7r^ as(/io) \ 512 vr y 

The expression (48) represents the "theoretical" side of the QCD sum rule. The "phe- 
nomenological" side is given by the r.h.s. of (47), with the (hadronic) spectral function Im 
ipis) written as 

ImV'(s) = lmip{s)\Res d{so - s) + lmip{s)\QCD 0{s - Sq) , (51) 

where the first term above describes the contributions of the resonances up to s = G.SGe^^ 
and the second term, i.e. the hadronic continuum, is identified as usual with the perturbative 
QCD expression, which in this case is given by 

-\rai,{s)\QCD = ^,mlis)s (l + 1^^) . (52) 

71 StT'^ \ 671 J 

The continuum threshold is expected to be close to the upper limit of the experimental 
data, i.e. sq — 6 — 7GeV^. In principle, though, sq is a free parameter. Predictions will be 
meaningful provided they do not depend strongly on the value of this parameter. 



Chiral dynamics provides a strong constraint on the behaviour of the hadronic spectral 
function near threshold, viz. 

^Im^(.) = ^MMI^a-^)!!-^), (53) 

with sf'^ = {Mk ± M^Y, and \d{s+)\ ~ 0.3 GeV^ A good fit to the experimental data 
1^ is obtained by using (53), which simulates non-resonant background, to normalize two 
Breit-Wigner forms for the -^"0(1430) and ii'Q(1950) resonances, with masses and widths: 
Ml = (1.40 ± 0.01) GeV, Ti = (325 ± 30) MeV, = (1.94 ± 0.03) GeV, T2 = (450 ± 100) 
MeV. 



As for the QCD parameters, we adopt the following values for the nonperturbative con- 
densates: ('Um)^o = —(0.25)^ GeV^ at a scale fio = 1 GeV and = 0.7 — 1. The 
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gluon condensate has been extracted some time ago from data on — e annihilation, 
and r-decay, with values in the range: (— G^) = 0.02 — 0.06 GeV^. The QCD scale for three 
flavours is A ~ 200-400 MeV 12^ ~ 



The invariant strange quark mass is determined by solving the equation resulting from 
inserting (51) on the r.h.s. of (47) and using (48) on the l.h.s.. Typical results are shown in 
Fig.3 (for A = 200 MeV), and Fig.4 (for A = 400 MeV), corresponding to {^C^)^^ = 0.02 
GeV^. In both these figures we used {ss)f^^J ('^'^)mo — 1- Results are essentially unchanged if 
this ratio deviates from unity by some 30%. The error on is determined by its variation 
when all relevant parameters are changed within the ranges indicated above. This gives, for 
the two extreme choices of A, 

rhs = 213 - 222MeV , m,(lGeV) = 171 - 179MeV (A = 200MeV) , (54) 

rhs = 142 - 147MeV , m,(lGeV) = 162 - 168MeV (A = 400MeV) , (55) 

where the variation of the strange-quark mass, for a given value of A, reflects the uncertain- 
ties in the values of the gluon condensate and sq. 



The results of this determination show a welcomed stability in the Borel variable M^, as 
well as in the continuum threshold Sq. To estimate the error induced by the uncertainties in 
the hadronic spectral function, we have varied the resonance parameters within the limits 
shown above. This gives an additional error of about ±7 MeV. The final uncertainty in 
is almost exclusively due to the influence of A and the gluon condensate. 

The effect of the three-loop radiative correction to IIo, and hence to has been to 
reduce the value of the invariant mass rhs by (5-10)%. Combining the results in (54) and 
(55) into a single prediction and including the additional error due to uncertainties in the 
hadronic parameters, leads to 

m, = 182 ± 45MeV , m,(lGeV) = 171 ± 15MeV. (56) 



4 Conclusions 

In this paper we have discussed in detail how to absorb mass singularities into the vacuum 
expectation value of other operators of appropriate dimension, for the case of the mass loga- 
rithms log(m^/(5^). We have also included the mixing with other dimension-4 operators to 
two- loop order. A comparison has been made with earlier analyses of this problem 0,[07 [0- 
0. In particular, we have shown that in our approach it is possible to avoid terms involving 
inverse powers of as which, being parametrically enhanced, might lead to large corrections. 
We have then used the QCD expression of the current correlator involving the strangeness 
changing vector current, together with a fit to the experimental data on the I = \, S-wave 
Kn amplitude, to determine the strange-quark mass through a Borel QCD sum rule. Our 
results for are in agreement, within errors, with the determination of which used the 
same fit to the data, but employed the QCD approach of to remove mass singularities. 
The errors we quote for the strange-quark mass are larger than those in 0. This is mostly 
due to the fact that in the gluon condensate was fixed at the single value (^G^) = 0.03 
GeV^, and A was allowed to change in the narrower interval A = 100 — 200MeV. 
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Figure Captions 



Figure 1. Lowest order contribution to the correlator iplQ'^). 

Figure 2. Vacuum diagram contributing to the pcrturbativc VEV of the operator ss. 

Figure 3. The invariant strange quark mass as a function of the Borel variable M^, for 
A = 200 MeV. The values of the gluon condensate (^G^)/io of the continuum 
threshold Sq have been varied between 0.02 and 0.06 GeV^ and respectively, Sq — 6 
and 7 GeV^. 

Figure 4. The same as Fig.3, except for A = 400 MeV. 

Figure 5. Results for the running strange quark mass 771.5(1 GeV) at the scale fi — 1 GeV. 

a) A = 200 MeV, b) A = 300 MeV, c) A = 400 MeV. The continuum threshold 
So = 6.5 GeV^. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409371vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409371vl 



Fig. 2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409371vl 
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This figure "figl-4.png" is available in "png" format from: 
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